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ABSTRACT 


Technical progress is reported on the work accomplished 
during the third quarter of contract No. DA 28-043 AMC-01253{E) 
entitled, Airborne Millimeter Radar Techniques Study (Obstacle 
Avoidance) with the United States Army Electronics Command, 

Fort Monmouth, New Jersey. During the reporting period, work 
continued on the following principal tasks: 

Advanced Sensor Techniques 

Mathematical System Modeling for Computer 
Simulation 

Millimeter Radar Experiments 
Advanced Sensor Techniques 

Research into clutter discrimination by means of fre¬ 
quency agility is presented. A model has been developed which 
exhibits the characteristics of observed clutter. The trans¬ 
mitted frequency has been varied and the electric field reradi¬ 
ated from the model, has been computed. The results indicate 
a marked dependence of the backscatter signal on the amount of 
frequency agility. Spectra are presented which indicate energy 
shift over the clutter band. 

The study of signal processing techniques for resolu¬ 
tion improvement has continued and preliminary results of an 
azimuth compression technique are presented. Optical data pro¬ 
cessing techniques have been investigated for use in identify¬ 
ing obstacles. 

Mathematical System Modeling for Computer Simulation 

A digital simulation of an intermediate frequency amp¬ 
lifier has been accomplished which incorporates limiting and 
bandpass filtering. A simple horn antenna has been simulated 
which can be used as the building block for more complex aper¬ 
ture illumination functions. 

Millimeter Radar Experiments 

The statistical analysis of radar data obtained at 
Ballistics Research Laboratory, Aberdeen, Maryland, has been 
completed and the probability density functions are presented. 

The data has been correlated with the power spectral densities 
published in ECOM-01253-2. 

The millimeter radar has been completed and is described 
in detail. Initial experiments have indicated that the system 
is very stable and can be used for reliable clutter and obstacle 
measurements. Data taken with the radar, at 70 GHz, is published 
and pertinent characteristics are discussed. 
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FOREWORD 


This program is being conducted under Contract Ho. 

DA 28-043 AMC-01253{E) with the United States Army Electronics 
Command, Fort Monmouth, New Jersey. The DA Task No. is 
1E1-20601-A-219-03-01. 

The work described in this third quarterly report 
represents significant progress toward achieving the goals of 
this program. Stated briefly, our purpose is to provide a 
program of study, investigation, and applied research which is 
intended to supply some of the basic information on which to 
base the design of future (1975 generation) airborne obstacle 
avoidance equipment for use in Army aircraft under tactical 
conditions. To meet these goal.*, our program is organized 
under three major technical tasks, namely: 

Task 20 - Physical Characteristics of Obstacles 

Task 30 - Advanced Sensor Techniques 

Task 40 - Mathematical Modeling for Computer Simulation 

Detailed descriptions of these tasks are given in Technical 
Report ECOM-01253-1. The remainder of this section describes 
the accomplishments of the work performed during this quarter, 
within the framework of the over-all pi >gram goals. 

Task 20 has been completed and has been described in 
detail in report ECOM-01253-1. Significant and useful data 
have been obtained on: 

a. surface density and height characteristics of 
vertical obstacles to be found in typical rural, 
suburban, and small urban areas 

b. attenuation and backscatter cross section charact¬ 
eristics of the atmosphere (including precipitation) 
at millimeter and optical wavelengths 

Task 20 also includes the following items, which are 
described in report ECOM-01253-2. 

a. an extensive annotated bibliography of references 
on the scattering characteristics of various types 
of target shapes 

b. an extensive tabulation of equations for calculat¬ 
ing the backscatter cross section of typical 
target shapes 

c. a study of possible target signatures 
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Task 30 continues the research into clutter discrimination 
techniques utilizing frequency agility. In order to evaluate 
the effects of frequency shifting techniques on the reradiation 
characteristics of scatterers, a suitable model was sought which 
displays the observed behavior of true clutter. Such a model 
has been developed, based on the multiple radiator model of 
Kerr [1] and Goldstein [2], and is shown to behave as actual 
clutter. The frequency of transmission is then varied in 
accordance with a sinusoidal function of time and the signal 
characteristics are observed. The results show a marked dependence 
of the backscatter signal on the amount of frequency agility. 
Spectra are computed which indicate the shift in energy over the 
band. 

Task 30 includes an initial investigation into optical 
data processing techniques which may be of use in identifying 
and recognizing obstacles. Such techniques have been developed 
for high resolution side-looking radar systems and may prove 
useful in the present study. 

The study of signal processing techniques for resolution 
improvement has continued and the work to date is described. 
Significant improvements in azimuth resolution may be achievable 
by the synthetic aperture approach. 

Under Task 40, a digital simulation of an intermediate 
frequency amplifier has been accomplished which incorporates 
limiting and band-pass filtering. A simple horn antenna has 
been simulated which can be used as the building block for more 
complex aperture illuminations and shapes. These programs are 
being combined with our previously determined clutter programs 
to develop a complete radar simulation. 

The statistical analysis of radar data obtained at the 
Ballistics Research Laboratory, Aberdeen Maryland, has been 
completed and the probability density functions are presented. 

The data has been correlated with the power spectral densities 
published in ECQM-01253-2 to furnish a complete statistical 
picture of several objects measured at 70 GHz. 

The millimeter radar has been completed and is described 
in detail. Initial experiments have indicated that the system 
is very stable and can be used for reliable clutter and obstacle 
measurements. Data taken with the radar, at 70 GHz, is published 
and pertinent characteristics are discussed. 
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TECHNICAL DATA 


1.1 Physical Characteristics oi Obstacles 

Task 20 has been completed. Reference should be made 
to reports ECOM-01253-1 and ECOM-01253-2 for the results of 
this investigation. 


1.2 Advanced Sensor Techniques 

Research into clutter discrimination by means of 
frequency agility is presented. A model has been developed 
which exhibits the characteristics of observed clutter. The 
transmitted frequency has been varied and the electric field 
reradiated from the model has been computed. The results in¬ 
dicate a marked dependence of the backscatter signal on the 
amount of frequency agility. Spectra are presented which 
indicate energy shift over the clutter band. 

Optical data processing techniques have been investi¬ 
gated for use in identifying obstacles. Further research into 
this area will be made in the next reporting period. The study 
of signal processing techniques for resolution improvement has 
continued and preliminary results of an azimuth compression 
technique are presented. 


1.2.1 Clutter Discrimination By Means F r equency Agility 


1.2.1.1 Introduction 


In order to evaluate the effects of frequency shifting 
techniques on the reradiation characteristics of scatterers, 
a suitable model was sought which displays the observed be¬ 
havior of true clutter. Such a model has been developed based 
on the multiple radiator model of Kerr [1] and Goldstein [2] 
and is shown to behave as true clutter. The frequency of trans¬ 
mission is then varied in accordance with a sinusoidal function 
of time and the signal characteristics are observed. The re¬ 
sults show a marked dependence of the backscatter signal on the 
amount of frequency agility. Spectra are computed which indicate 
the shift in energy over the band. 


1.2.1.2 The Scatterer Model 


The development of a model to obtain the desired power 
spectrum characteristics was ach.eved through the use of a group 
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of randomly moving scatterers. As shown in Figure 1, "n" 
scatterers are arranged such that their mean positions are 
uniformly distributed in range. Each scatterer is permitted to 
move at an arbitrary velocity and to have an arbitrary starting 
position. The resulting time function observed at the receiver 
is the vector sum of the electric field from each of the scatter¬ 
ers. In the sections which follow, the results of such a model 
will be examined and compared with predicted clutter characteris¬ 
tics. 


1.2.1.3 Clutter Characteristics 

The signal radiated by a large group of scatterers has 
been studied on a statistical basis for many years. The 
theoretical models which have served to describe observed pheno¬ 
mena generally consist of an infinite number of independent 
scatterers, each free to move at a different velocity. The 
description of the total signal emanating from such an array 
forms a branch of the general theory of random processes dating 
back to Rayleigh [3]. 

The total received electric field strength can be shown 
to be of the form [4]: 


E(t) - R(t)cos[u> o t + 


(t) ] 


( 1 ) 


where is the transmitted frequency and R(t) and <f(t) are 

functions of scatterer motion. If we define a vector $(t), whose 
magnitude is equal to R and whose phase is 4>, then the probabil¬ 
ity density function of R, i.e., the probability that R lies 
between R and R + dR, is given by the well known Rayleigh 
distribution [5]: 

R 


p bF) ■ 2 (-J=) 


-e=y 


o o 

A plot of equation (2)is given in Figure 2, 


( 2 ) 


The mean value 


/r 


of-is 

/p~ 


and the rms value is 


/ ( 4“ 71 ) . 


2 If we consider the signal energy, P, which is proportional 
to R , then the probability density function of P is given as [5]: 
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A plot of equation (3) is shown in Figure 3. 


Extending the discussion to the case of a fixed target 
return plus clutter, Schwartz [6] has shown that, 

-£-♦«*> 


- (■ 


2Re 

p(R) » p— 
o 


I o (2mR//P^) 


(4) 


where I is the modified Bessel function, and m 2 is the signal- 
to-noise power ratio. For m 2 >>l, equation (4), reduces to the 
Gaussian function: 


p(R) 



(5) 


Thus, for high signal-to-noise ratios, the power spectral^ density 

rapidly approaches a Gaus sian distribution with mean m/P and 

standard deviation /P /2. ° 

o 

Up to this point we have been concerned with the 
long-term statistics of clutter which essentially determine 
the averages [6] of the noise process. Now let us examine the 
correlation between samples of the clutter signal spaced at t 
second intervals, i.e., the autocorrelation function of the 
clutter process. This problem has not been solved analytically 
for the model discussed previously although much experimental 
data has been obtained for a wide variety of scatterer mobilities 
and forcing functions, such as wind. The Fourier transform of 
the autocorrelation function is the power spectral density and 
it is with this function that we shall be concerned. Barlow [7] 
suggests that a good approximation to the experimentally 
determined spectral density, valid over a wide range of radar 
frequencies, is given by: 


W(f) 



( 6 ) 


where "a" is a spectrum parameter dependent on the composition 
of the area being observed and f is the frequency of operation. 

In terms of these two quantities, namely the probability 
density function and the power spectral density, we can now 
discuss the random scatterer model. 
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1.2.1.4 Experimental (Computer) Results 


The scatterer model discussed In Section 1.2.1.2 was pro¬ 
grammed for n»10 and has been run at both 16.5 GHz and 70 GHz 
without frequency agility. The resultant time function is shown 
in Figure 4 Detailed analyses of the nature of the time function 
were made to determine its reliability as a realistic clutter 
simulator. 

Figure 5 shows the power spectral density of the model 
while Figure 6 shows the probability density function of the 
model, both taken at 16.5 GHz. It is seen that the statistics 
closely resemble those obtained by Barlow and others [7] for 
actual terrain. In addition, since the velocity of the scatterers 
is arbitrary, any width spectrum can be obtained merely by 
scaling the time base or conversely, the frequency base. 

Figures 7,9, and 11 show the resulting time functions 
when frequency agility is applied to the model. The statistics 
of these functions are plotted in Figures 5,8,10 and 12. 

The similarity between the signal statistics from the 
model and those obtained experimentally is quite good if 
reference is made to Kerr [8]. Figures 13 and 14 show the 
measured probability densities of wooded terrain taken at 
10 and 25 mph wind velocities. It should be noted that the 
probability density goes from exponential to Rayleigh as the 
wind velocity increases. Comparison between these curves and 
Figures 6 and 10 shows that frequency agility on the present 
model produces the same effect as increased wind velocity. It 
can be seen that an abrupt change in clutter power, or standard 
deviation, occurs as the magnitude of the frequency shift in¬ 
creases beyond 4 MHz. In addition, pockets of energy appear 
around harmonics of the shift frequency. A close examination 
of the time waveform indicates that the high frequency energy 
is a result of intense interference and occurs infrequently 
even with frequency agility. The results of these experiments 
indicate that very large amounts of frequency agility achieve 
little in whitening the clutter spectrum beyond a certain point. 
This critical frequency is quite abrupt and will be examined 
closely in subsequent studies. In addition, slowly moving 
objects will be placed in the clutter and their characteristics 
examined for possible signatures. 


Experimental data at 70 GHz will be presented in the final 
report. 
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. Magnitude of the Square of the Electric 
Field Reradiated by a 10-Scatterer Model. 
Frequency of Operation, f, is 16.5 GHz with 
a Pulse Repetition Frequency, PRF, of 2400 
Hz and Zero Frequency Agility. The Separa¬ 
tion between Scatterers is 20 Feet. 
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Power Spectral Density, W(f), of the 
n-Scatterer Model, for n = 10. Af is 
the magnitude of frequency shift. 
























6. Probability Density of the Signal of 
Figure 4. N is the Number of Samples 
Recorded. 
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Figure 7. Magnitude of the Square of the Electric Field 
Reradiated by a 10-Scatterer Model. Frequency 
of Operation is 16.5 GHz with a Pulse Rep¬ 
etition Frequency, PRF, of 2400 Hz. Frequency 
Agility with a Magnitude of 4 MHz and a 
Slew Rate of 60 Hz was Employed. The Sep¬ 
aration between Scatterers is 20 Feet. 
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8. Probability Density of the Signal of 
Figure 7. N is the Number of Sample 
Recorded. 
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Figure 9. Magnitude of the Square of the Electric Field 
Reradiated by a 10-Scatterer Model. Frequency 
of Operation is 16.5 GHz with a Pulse Rep¬ 
etition Frequency, PRF, of 2400 Hz. Frequency 
Agility with a Magnitude of 10 MHz and a 
Slew Rate of 60 Hz was Employed. The Sep¬ 
aration between Scatterers is 20 Feet. 
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Figure 11. Magnitude of the Square of the Electric 

Field Reradiated by a 10-Scatterer Model. 
Frequency of Operation is 16.5 GHz with a 
Pulse Repetition Frequency, PRF, of 2400 Hz. 
Frequency Agility with a Magnitude of 20 MHz 
and a Slew Rate of 60 Hz was Employed. The 
Separation between Scatterers is 20 Feet. 
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Figure 14. 



Probability Density Function of Ground 
Clutter Taken at 9.2 Cm. Heavily Wooded 
Terrain, Wind Speed 25 MPH. P is the 
average power. {8] 












1.2.2 Resolution Improvement By Means of Signal Processing 




1.2.2.1 Introduction 

High resolution in both range and azimuth angle can be 
achieved simultaneously by means of signal processing techniques. 
High resolution is achieved utilizing existing low resolution 
equipment/ by the addition of a frequency modulation capability 
to the transmitter and an appropriate signal processing capabil¬ 
ity to the receiver. Although a long pulse and a wide azimuth 
beam are used, the range resolution of a short pulse and the 
azimuth angle resolution of a narrow azimuth beam are realized. 

An azimuth compression technique is used analogous to 
the commonly employed pulse compression techniques. Frequency 
modulation is introduced under the azimuth space pattern. 

The carrier frequercy of each individual pulse is kept constant, 
but the carrier frequency is varied from pulse to pulse. By 
means of suitable processing, a compressed azimuth space pattern 
is realized. 

Pulse compression techniques commonly transmit a pulse 
envelope with a linearly increasing carrier frequency. The 
received pulse is applied to a filter that has a delay charact¬ 
eristic which is a linearly decreasing function of frequency. 
Thus, the lowest frequency which is received first is delayed 
the longest, and the highest frequency which is received last, 
i* delayed the least, resulting in a compressed pulse. 

An analogous technique is employed in this report for 
compression in azimuth. Varying the carrier frequency from 
pulse to pulse introduces frequency modulation within the 
azimuth space pattern envelope in a manner analogous to the 
introduction of frequency modulation within a pulse envelope. 

The series of pulses of different carrier frequency are applied 
to a filter, with delay varying with frequency. The filter de¬ 
lays each pulse the appropriate amount so that the pulses are 
all added at the same time, resulting in a compressed azimuth 
space pattern. 

However, there are important differences between using 
a pulse with modulated carrier for range compression and using 
a group of pulses with varying carriers for azimuth compression. 

A pulse with a linearly modulated carrier, as it travels 
out in range, illuminates, in turn, each of two point targets 
separated in range. The returns from both targets are pulses 
with carriers linearly increasing in an identical manner. 

However, the two targets do not return the same frequency at 
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any given time, even though they are separated by less than a 
pulse width. Tims, each frequency is returned from the two 
targets at a different time, and one compression filter is 
sufficient to separate the returns. Since the compression 
process is a linear one, the output consists of two compressed 
pulses. 


The situation is quite different for a series of pulses 
with varying carriers. As the antenna scans in azimuth, a 
target may pass into the beam at any time. Although the return 
from a target consists of a series of pulses with one complete 
cycle of the carrier variation* 2 } the sequence in which the 
carrier frequencies are returned depends upon where in the cycle 
the return begins. Thus, the return from a target is a series 
of pulses which are frequency coded; the code being the sequence 
in which the carrier frequencies are returned. 

One compression filter, matched to that particular 
sequence, is required for each sequence. The filter for each 
sequence has an output only when it receives the correct 
sequence. The filter delays each pulse of the sequence an 
appropriate amount and then adds them at the correct time. 

When two point targets are in the azimuth beam, they 
both return the same frequency. Two point targets will only 
return differing sequences if they are separated in azimuth by 
at least the angle that the antenna moves in the time between 
pulses. This is the azimuth resolution limit of this technique. 

The summation of all the pulses of a sequence at the same 
time also achieves pulse compression. The summation is that of 
a number of pulses with varying carriers. For a linear variation, 
the result is a compressed pulse of width equal to the reciprocal 
product of the number of pulses and the frequency variation from 
pulse to pulse. 

High azimuth resolution may also be achieved by doppler 
techniques. However, a fundamental shortcoming of the doppler 
method is that resolution is a function of the angular position 
of the targets. This high resolution technique fails for targets 
on or near dead ahead. The technique of varying the carrier fre¬ 
quency from pulse to pulse does not have this shortcoming. 


This simple case is used for purposes of illustration. 
Modifications are necessary in a practical situation. 
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Section 1.2.2.2 gives an explanation of the mathematical 
foundation upon which the range-azimuth compression technique 
is based. Section 1.2.2.3 gives an explanation of how high 
range and azimuth angle resolution are realized, in particular, 
by means of a varying carrier from pulse to pulse, based upon 
the mathematical foundation of section 1.2.2.2. 


1.2.2.2 Combined Range-Azimuth Angle Ambiguity Function 

In order to be able to resolve two targets, the signals 
returned from the two targets must be separated on the basis of 
their differences. It an be shown that a measure of the 
differences between the signals returned from two targets, one 
at range tj,( range given in terms of the time delay of return, t), 
and azimuth angle 8j, the other at range t 2 and azimuth angle 
62 is the combined range-azimuth angle ambiguity function given 
by (9] 


J t,v ( t ' v) " J t (t)J v (u o v) 


(7) 


where v (t,v) » combined range-azimuth angle ambiguity function 

J (t) ® range ambiguity function ■ autocorrelation 

1 function of pulse * /a(t) a*(t-x) dt, where (8) 

star denotes complex conjugate. 


J v U 0 v) ■ azimuth angle ambiguity function « autocorrelation 
function of azimuth space pattern = 

I G*(wqW)G t“ 0 (y + v)]dy (9) 

s(t) « returned signal (3) 

= a(t) e^“o t 


a (t) 


03 


o 


T 

G(w 0 p) 

v 


u + v 


= shape of pulse 
® carrier frequency 

= T 2 " T 1 

® azimuth space pattern as a function of angle 
sin 6i 
- 

sin e 2 
c 
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A signal s(t) may be represented in 
s(t) = a(t)e jlui o t " * (t)J = a(t) 
where a(t) = a(t) e^*^ 


complex notation as 

ioj t 

e o 
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It can be seen from equation (7) that the combined 
range-azimuth angle ambiguity function is separable into the 
range ambiguity function and the azimuth angle ambiguity 
function. Thus, range resolution and azimuth resolution are 
independent of one another. 

The range ambiguity function is the autocorrelation 
function of the pulse and the azimuth angle ambiguity function 
is the autocorrelation function of the azimuth space pattern. 
The range resolution properties are given by the shape of the 
pulse autocorrelation function and the azimuth angle resolution 
properties by the shape of the azimuth space pattern autocorre¬ 
lation function. High resolution is achieved if signals with 
suitable autocorrelation functions can be found. It will be 
shown in section 1.2.2.3 that the signal, with carrier varying 
from pulse to pulse, has an autocorrelation function that gives 
.both azimuth space pattern compression and pulse compression. 

The output of a filter matched to a signal is the 
autocorrelation function of that signal as a function of time. 
The output of each filter matched to a sequence of pulses is 
the autocorrelation function of the sequence to which it is 
matched. 


1.2.2.3 Compression in Range and Azimuth by Means of Frequency 
Coding from Pulse to Pulse 

Based upon the mathematical foundation of section 1.2.2.2, 
an explanation is given of range and azimuth angle compression 
by means of frequency coding. 


4 If s(t) is a time functii whose Fourier Transform S(w) = 

a(oo)e-^ B , the response of a filter, h(t) , matched to this 
signal is given by 

h(t) = s*(-t), where * denotes complex conjugate 
H(u>) = S*U) = a(u))e~ jB(u,) 

The output of the matched filter, e(t), is given by 

e (t) = Jv / S(u)S*(aj) e^ ut dw « i / a 2 ( u )e jll,t d w 

ZTT J Z IT ] 

By adding all the frequency components in phase, the matched 
filter produces a compressed pulse at its output. 


* 
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1.2.2.3.1 Azimuth Compression 


The azimuth space pattern is converted to a function 
of time at the receiver by means of antenna rotation to give 


G[u> o ii (t-y) 1 


where 


= carrier frequency 
u = antenna scan rate 

y = time at which target passes into beam 

A typical sequence of pulses returned from a target, 
with carriers varying linearly from pulse-to-pulse, is given 
by (6 » 5 6 7 > 

s (t ) = ^ 1+m gu &t)e jl “o +(n ~ m,u i 1 lt-(n-m)T-t ] 

o L * o 

n=m 

t-(n-m)T-t 

rect ( ---) (11) 

P 

where s R (t) “ sequence of pulses, n * 1,2,3,... 

m = pulse at which sequence begins 
T = 1/PRF 
ip = pulse width 
T t = delay of target in range 
uj = carrier change from pulse-:o-pulse 


rect (z) = 1, 


0< z -1 


5 The effect of the variation of carrier frequency from pulse- 
to-pulse on the azimuth space pattern is neglected, since 
the variation is small compared to the original carrier. 

6 The phase of the azimuth space pattern is neglected. g(u) 0 ut) 
represents the amplitude of the azimuth space pattern as 

a function of time. 

7 For purposes of illustration, it is assumed that the carrier 
is linearly varied and that the return from a point target 
contains only one cycle of the carrier variation. 
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The azimuth space pattern is considered constant during a pulse 
time giving 

N-l+m t-(n-m)T-t 

s (t) = l glw vi(n-m)T] rect (----) 

n=m P 

( 12 ) 

e j U o +(n-m)w 1 l (t- {n-m)T-t t ] 

The pulses of the sequence are delayed appropriately 
and then summed at the same time by the filter matched to that 
sequence. The output of the matched filter, m(t), is given by: 

m(t) = e^o*' ^ m g U p (n-m)T) rect(^-) t (n-m) uj) t’ (13) 
n=m T p 


where 


t' = t-(n-m)T--r t 


1.2.2.3.2 Range Compression 

If the pulses are all of equal amplitude, the filter 
output given by equation (13) becomes: 


m' (t) 


= e 3- 0 t' N £ 1+m e j ((n-m)u,]t' 


n=m 


rect (j -) 
P 


(14) 


N . , 
sin — I.) 1 1 


= l “o + u ’i J t ' 


“i 

sin y t' 


rect(^—) 
‘ P 


Equation (14) is a frequency shifted carrier with envelope 


sin 


sin 


N 

2 


u x t' 



rect(^—) 
T P 


(15) 
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Equation (15) is a compressed pulse of height N and width -jj 
T husi pulse compression is realized. The results for 
N-10, t » .2 nsec, f=4MHz, f =100MHz are given in Figure 15. 
P ^ 

1.2.2.3.3 Autocorrelation Function 


It can be shown that equation (13) is the center term 
of the autocorrelation function (the other terms are eliminated 
by the filtering process) of a train of rectangular pulses with 
carriers varying linearly from pulse to pulse [10). Thus a 
filter with this output is a close approximation to a matched 
filter for that particular frequency-coded sequence. 


1.2.2.3.4 Unequal Amplitudes and Random Phases From Pulse to Pulse 

Equation (13) is the sum of N terms of a Fourier series. 
Equation (14) gives the form of the sum when all the terms have 
equal amplitudes and equal phases. It is desired to determine 
the form of the sum when the amplitudes and the phases vary 
from term to term. 

If the terms have unequal amplitudes, but equal phases, 
the summation is that of N terms of a Fourier series with real 
coefficients. The terms all add in phase. The coefficients 
are given by the azimuth space pattern sampled at regular in¬ 
tervals. If the amplitudes are not greatly different from one 
another, the envelope is still of the form of equation (15) given 

by: N 

sin tt Wl t' 

k --- (16) 

“■ 1 

sin —- t' 

However, k is smaller than 1. An example is given in Figure 16 
which is identical to Figure 15, except that the pulses have 
random amplitudes. However, if the amplitudes of certain 
components predominate, these determine the shape of the en¬ 
velope. 


The form of the sum when the phases vary determines 
the degree of coherence required in the system. If the terms 
have equal amplitudes, but the phases vary, the summation is 
that of N terms of a Fourier series with imaginary coefficients. 
In the worst case, with the phases varying in a random manner, 
the terms do not add in phase and the form of the envelope of 
equation (16) is not retained. Figures 17, 18, and 19 show the 
summation of terms of random phases. The random phases of 
Figure 17 are multiplied by 45° and 90°, respectively, to obtain 
the phases of Figures 18 and 19. Appreciable degradation of the 
envelope does not appear until the values of the random phases 
approach 45°. 
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Figure 15. Sum of 10 Pulses with Linearly Varying 
Carrier from Pulse to Pulse. The Sum 


9 

2 sin [2w(f + nf.)t] is Evaluated with 

n=0 ° L 

f = 100 MHz, f. = 4MHz and r = 0.2 us. 
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Figure 17. Sum of 10 Pulses with Linearly Varying 

Carrier of Random Phase from Pulse i'u' se. 
The Random Phases <f> n are Less Than 1". Tae 
Sum 9 

£ cos [2 it (f + nf^)t + $ n ] is Evaluated 
n=0 

with f = 100 MHz, f. = 4 MHz and t = 0.2 us. 
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Figure 18. Sum of 10 Pulses with Linearly Varying 

Carrier of Random Phase from Pulse to Pulse 
The Random Phases $ n are Less than 45°. 

The Sum 9 

Z cos [2 n (f + nf,)t + $ n ) is 
n=0 

Evaluated with f = 100 MHz, f^ - 
and = 0.2 us. 


4 MHz 
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Figure 19. Sum of 10 Pulses with Linearly Varying 

Carrier of Random Phase from Pulse to Pulse 
The Random Phases 4 are Less than 90°. 

The Sum ^ 

E cos [2ir(f + nf.)t + $ ] is 
n =0 o l n 

Evaluated with f = 100 MHz, f. = 

C X 

and t = 0.2 ys. 
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1.2.3 Electro-Optical Information Processing 

The high data rate required for low-altitude obstacle 
avoidance necessitates a close examination of data processing 
and display techniques which can handle large amounts in real 
time. In this section we consider the basic principles of 
electro-optical data processing which is playing an increasing 
role in high resolution three-dimensional sensors and displays. 
The theory underlying optical processing has been studied for 
a number of years and is being examined for its ability to 
process large amounts by information for use in obstacle avoid¬ 
ance applications. 

1.2.3.1 Introduction 


Because certain useful filtering operations may be 
carried out with greater facility by optical means than by 
electronic means, the application of electro-optical techniques 
to the processing of information in a radar system appears to 
offer certain advantages. One of the primary advantages of 
optical information processing systems is the high information 
handling capability resulting from two-dimensional processing 
as opposed to the one-dimensional processing possible in an 
electronic circuit. In addition, large quantities of information 
may be processed efficiently because one process may be simul¬ 
taneously applied to many signals. Similarly, several processes 
may be simultaneously applied to the same signal. This results 
in reduction of system size and complexity. 


1.2.3.2 Optical Data Processing 

Two general groups of optical data processing are or 
interest, coherent and incoherent. In both groups, input-output 
relationships analogous to electronic processors can be obtained. 

The first group, the coherent processors, operates on 
the signal to form its Fourier transform optically. This 
Fourier transform, which is actually a diffraction pattern, can 
then be operated upon by spatial techniques. Thus, the coherent 
processors are equivalent to electronic filters and consequently 
any linear data processing operation can be performed in a 
coherent optical processor. 

The basic space-coherent processor, shown in Figure 20, 
consists of two lenses with their optical axes and focal points 
coincident, illuminated by a space-coherent light source such 
as a laser. (Since only spatial coherence is required, a mono¬ 
chromatic pin-hole source and suitable collimating equipment 
could be used.) The light source impinges upon the "signal" 
which is located at the object plane of the first lens. This 
signal can be a photographic transparency, or in the case of 
real time systems such as radar, an ultrasonic light modulator. 
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Figure 20. Basic Space Coherent Information 
Processor {II) 
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An ultrasonic wave passing through the modulator causes small 
variations in the index of refraction of the material. As a 
result, the phase of the entrant plane wave of light will be 
retarded in some portions of the modulator and advanced in 
others. The plane wave then emerges as a modulated wavefront. 
It is expected that these modulators will eventually be fab¬ 
ricated with frequency responses of several gigahertz. The 
first lens then forms the diffraction pattern of the signal in 
the diffraction plane between the two lenses. The light amp¬ 
litude and phase of the diffraction pattern are the Fourier 
transform of the signal in the object plane. If the amplitude 
function is then squared (analagous to an electrical power 
spectrum), the light intensity as a function of position is 
obtained. 


Processing of the signal is accomplished by simple 
spatial filtering (masking) of the diffraction pattern. The 
light then passes through the second lens which performs 
another transform and the processed information is displayed 
in the observation plane of the second lens. 

The second group of optical processors, the incoherent 
processors, have a somewhat more simple configuration than the 
coherent processors. A basic incoherent processor, shown in 
Figure 21, consists of a large diffuse light source, two signal 
sources (or a signal and a reference source) and a lens, arranged 
in that order, with their axes coincident. The signal sources, 
which define two object planes for the lens, can be ultrasonic 
light modulators or simply photographic transparencies. As 
rays of light leave the diffuse source at some angle, e, with the 
optic axis, they pass through the first modulator (object plane 
1 ), and this light acquires the distribution of the first 
signal fi(x,y). The rays then pass through the second modulator 
(object plane 2) and the distribution of the light is fj(x,y)- 
fp (x + x q , y + y ), where f ? (x,y) is the second signal and x q 
and y Q are functions of b. The lens then integrates this 
light and, in .the observation plane of the lens, the distribution 
is given by: 



fi(x,y)f 2 (x+x Q ,y+y Q )dxdy 


(17) 


This is seen to be a two dimensional analog of the finite time 
cross-correlation function 111]: 


J fi(t)f 2 (t+r)dt 
-T 


( 18 ) 
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Figure 21. Basic Space Incoherent Information 
Processor [11] 


34 


9 

















niittifriiiHiiinwHiiwnwMli 1 






It should be noted that in the system just described, 
the diffuse light source and the first signal source may be 
replaced by a cathode ray tube. 

Because signals are processed in two dimensions in the 
incoherent processor (as was the case for coherent optical 
processors ), the information packing density is quite high. 

Many signals may be processed in parallel without the size or 
complexity of the system increasing. 

One feature of the optical correlator which is of 
interest in radar range measurements is space delay. If the 
second signal f 2 (x,y) is shifted in the x or the y plane, a 
proportional shift will be observed in the correlation function 
at the observation plane. The delayed signal from the radar is 
recorded as a proportional space delay in f 2 (x,y). The position 
of the correlation function peak determines the time delay. 


1.3 Mathematical System Modeling for Computer Simulation 


1.3.1 Intermediate Frequency Amplifier Simulation 


1.3.1.1 Introduction 


This task, is a continuing effort to obtain a digital 
computer simulation of a radar system. This section discusses 
the simulation of a linear IF amplifier in cascade with several 
nonlinearities. 

1.3.1.2 Mathematical Development 

The IF amplifier is shown schematically in Figure 22. 

The unity gain amplifiers are shown to indicate impedance 
isolation. The transfer function of the tank circuit is written 
down from inspection: 


v i RLCs" + sL + R 
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The parameters R, L, C are computed using equations 



= 1/LC 

(20a) 

Q 

“ R/eL 

(20b) 

R 

= 1000 ohms 

(200 


Q and <u are designer input quantities and the value of 
R » 1000 was arbitrarily chosen as a practical value for 
transistor circuits. Equation (20a) is an approximation and 
may be in error by about 1% for low-Q circuits. This was 
compensated for in the program. 

For the steady state frequency sweep, vector algebra is 

used to calculate the magnitude and phase of the first stage output 

for a unit amplitude sinusoidal input. This output is used as the 

input for the next stage. No nonlinearities are included and 

the input freouency is swept from . 5f__ to 1.5f T _, where 

IF IF 

f IF = the intermediate frequency. The results are plotted 

by computer and are shown in Figure 23. No attempt was made 
to optimize the parameters for any specific application. 

The input to the amplifier for the transient analysis 
was a sinusoid plus noise, such that, 

input = 0, t<0 

input = sin a> in t + noise, taO (21) 

The input frequency is a read-in quantity and the noise 
is band unlimited with a flat distribution. The problem is 
allowed to run for 10 cycles of the input frequency. The 
differential equations for the filters are: 

RLC v* = Lv. - Lv - Rv 
o i o o 


v = f v dt 
o Jo 



There is no unique integration procedure fox solving 
these equations on a digital computer. Based on past experi¬ 
ence, the input sinusoid was divided into 60 parts (.5 nanosec. 
in width for a 30 MHz IF) and a combined Euler and Modified 
Euler integration algorithm was used. Equations (21) and (22) 
were evaluated and the results stored on a magnetic disk. A 
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FREQUENCY IN MHZ 


Figure 23. Steady State Frequency Response of a Two 
Stage IF Amplifier Without Limiting for 
fj f =30 MHz; f x = 35 MHz; Q x = 5; 

f 2 = 25 MHz; Q 2 = 6 
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logic subroutine limited the input to the second filter. 

Figure 24 is a plot of the noisy input. Figure 25 is the output 
of the first filter, Figure 26 is the output of the second filter. 
All the graphs are plotted to the same scale and can be super¬ 
imposed. After the computation is complete, the stored values 
are _atrieved from the disk. The program stores the maximum 
and minimum values of the outputs and automatically scales the 
plot. 

1.2.2 Rectangular Aperture Simulation 


1.3.2.1 Introduction 


This section describes the digital simulation of the 
radiation pattern emanating from a rectangular aperture. The 
aperture can be a waveguide feedhorn or a linear planar array. 
The results of this simulation can be used as a design tool 
and also serve as the first step in the analysis of a multi¬ 
element (lenses and reflectors) antenna. 


1.3.2.2 Mathematical Development 


The field strength at a point 
radiating source is given by [12], 

_ lii 


as a result of a 


Up = 


f 


FU.n) 


<¥ + ?>V r l + T 2 V* 


d£d 


(23) 


See Figure 27 for the definition of coordinate systems and 
variables. If the far field pattern is of primary interest, 
equation (23) can be simplified to: 


U = 
P 


j2TTR 

X 


XR 


/ 


F { C , n) e 


2 7T ' 
X 


sin9(£cos<t> + rising) 


d£dn 


(24) 


Physically, equation (24) assumes that: 



b. All regions of the aperture are excited in time phase 

c. Most of the energy is concentrated near the z axis 
where cos 8 = 1 
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TIME IN MICROSECONDS 


Noisy Sinusoidal Input V = sin wt + noise 
for White Noise of Zero Mean and Unity 
Variance 
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Figure 25. Output of First Stage of IF Amplifier 
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AMPLITUDE IN VOLTS 
















Finally, if only the relative pattern amplitude is of interest, 
the equation reduces to: 


U <e,$) 


b/2 

■/ 


-b/2 



FU, n)e 



sine(£cos$ 


+ 


nsin$) 

d£dn (25) 


Figure 28 shows the spherical coordinate system used. e, and <f> 
are normal spherical coordinate angles with $ representing a 
rotation about the normal to the aperture plane and 8 a rotation 
with respect to 2 about the Y' axis. When 4 - 0°, 8 is in the 
principal plane, XZ (H plane), and when $ = 90°, 8 is in the YZ 
plane, (E plane). 


1.3.2.3 Computational Procedure 


Us 1 ng some relatively simple expressions for F(f.,n), 
equation (25) can be integrated and solved in closed form. It 
is these solutions which served as a check on the computer 
evaluation. Of course, a computer performs quantized numerical 
integration but is able to evaluate almost any complex integrand 
and sum the integrand with vanishingly small error. When the 
computer procedure cnecks a closed form integral evaluation, 
it can reasonably be inferred that the computer program is valid 
for integrands which cannot be solved in closed form. 

The procedure used in this problem was to conceptually 
visualize the aperture partitioned into small squares (see 
Figure 29). At the center of each square, a one-way isotropic 
radiator was placed and scaled so that the total energy and 
energy distribution in the aperture was unchanged. After an 
angle $ was selected, the integral was evaluated at 50 values 
of 6 (from 0 e - 5a/x). 


1.3.2.4 Evaluation of Results 


Figures 30 and 31 are computer plotted field patterns 
of a rectangular aperture with a/ a = .1. The curves are 
identical except that Figure 30 has a linear voltage ordinate 
and Figure 31 is plotted in db. The aperture excitation is 


F(£, n) = cos i) 


(26) 
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ANGLE IN RADIANS 


Figure 30. 



H-Plane Field Intensity Pattern. 
Cosine Excitation with = 0.1, and 

d 

Linear Scaling. 






























which is the easiest way to excite an aperture. The excitations 
are uniform (independent) with n. All output values are 
normalized to unity or zero db. The following table shows 
the excellent agreement between the results achieved and the 
textbook values. 


Table I 


Computer Results 


Textbook Values 



Half Power 
Beam width 



P'. ase distort, ion 


.06 rad 
(.6 A/a) 

.6 A/a 

.15 rad = 1.5 A/a 

1.5 A/a 

• 23 db 

-23 db 

(zero, odd lobes 
(f, even lobes 

(zero, odd lobes 
jtr, even lobes 

(from printed 
listing) 


.105 rad 
=1.05 x/a 

A/a 



Separation 
o f zeros 


Figures 32 through 35 are representative of a small 
rectangular waveguide feeding a large dish. The following 
information can be extracted from these curves. 

a. In the "H" plane, the beam width is 1.2 A/a which 
agrees with standard textbooks. 

b. In the "E" plane, the beam width is .94 A/a 
(correct value = .88 A/a). 

c. The first zero in both planes occurs at 90°. 

d. However, the energv distribution av large values 
of 0 appears to ha incorrect. The reason for this 
error is that the original equation assumed that 
the majority of the energy was concentrated normal 
to the aperture. This assumption is not valid for 
small apertures. 
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ANGLE IN RADIANS 


Figure 32. Voltage H-Plane Pattern for a Cosine 
Excitation with —= 1, and Linear 

Cl 

Scaling. 
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ANGLE !N RADIANS 


H-Plane i'uitenna Pattern for $ 
Cosine Excitation with —— •= 

cl 

db Scaling. 














0 .2 .4 .6 .8 1.0 1.2 

ANGLE IN RADIANS 


Figure 34. E-Plane Antenna Pattern for $ = ir/2. 

Uniform Excitation with —= 1, and 
Linear Scaling. 
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Figure 35. E-Plane Antenna Pattern for $ = n 
Uniform Excitation with --g— = 1, 
db Scaling. 
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2 . 3.2.5 S uiw.a r v 


T!'c concept of computing the radiation pattern of a 
radiatinc; attj by subdividing it into a large number of equiv¬ 
alent isotropic radiators appears to be valid. This simulation 
should be considered as only the first step and can be extended 
to include small apertures, circular apertures, phased arrays, 
reflectors, lenses, and multiple elements. 


1.4 Expe rimental Millimeter Radar 

The statistical analysis of radar data obtained at 
Ballistics Research Laboratory, Aberdeen, Maryland, has been 
completed and the probability density functions are presented. 

The data has been correlated with the power spectral densities 
published in ECOM-01253-2 to furnish a complete statistical 
picture of several objects at 70 GHz. The millimeter radar has 
been completed and is described in detail. Initial experiments 
have indicated that the system is very stable and can be used 
for reliable clutter and obstacle measurements. Data taken 
with the radar at 70 GHz is published and pertinent chiracteris- 
tics are discussed. 


1.4.1 Radar Design 

Construction of an active radar system, operating at 
70 GHz was completed during the past quarter. The system has 
undergone tests and is currently being used for measurements 
being performed in support of the study program. Portions of 
the system and the progress in its construction were described 
in Technical Reports FCOM-Ot?53-1 and ECOM-01253-2. Some svstem 
design details not described before, as well as actual system 
test data, are reported here. A photograph of the system is 
given in Figure 36. 

Significant improvements in the fabrication techniques 
for the BL-234C magnetron were evolved by the design engineers 
at the Pomac Division of Varian Associates during the past 
quarter. Failure analyses performed by Bomac on earlier tubes 
delivered to Harry Diamond Laboratories and to Norden, have led 
to revised cathode and magnet activation techniques. The magnetron 
finally delivered for use in the Norden system has exhibited 
excellent characteristics m te^mr cf mode width and stability. 

The peak power output of the magnetron was measured in excess 
of 900 watts. The frequency of operation was measured to be 
70.42 GHz. 
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The BL-S-517 crystal-protecting TR tube, also manu¬ 
factured by Bomac, was tested prior to installation in the 
system and the following significant data was obtained: 

Recovery time: approx. 0.2 microsecond 

Spike energy leakage: 0.04 erg 

Insertion loss on receive: 1 db 

"ith no apparent crystal degradation over the first few hours 
of system operation, it appears that the protection offered 
by the BL-S-517, is adequate for the particular 1N2792 crystals 
currently in use. 

Waveguide losses in the transmit direction from the 
magnetron flange to the input of the antenna feedhorn run were 
found to be 3.1 db. The losses in the receive direction from 
the feedhorn flange to the mixer input were measured at 4.2 db. 
The feedhorn run adds approximately 0.6 db to the system losses 
in each direction. 

Some experimentation with a simple choke design will be 
carried out in an effort to reduce the system losses at each 
flange junction. Nine such junctions are included in each 
siqnal direction. 

The system contains a variable PRF generator with 
outputs at 2, 3, and 5 kHz. An externally controlled PRF 
option is also available. Solid state devices are used through¬ 
out the PRF generator as well as in the system internal power 
suoplies. Relay provisions for warm-up time delay, modulator 
protection, ana magnetron filament power reduction in the trans¬ 
mit mode have been made. 

A relatively novel approach was taken in the design of 
the system modulator. Some experiments were performed with a 
free-running blocking oscillator which pulsed the magnetron 
directly. Tests revealed this scheme to be lacking in flexi¬ 
bility and reliability, particularly when employed with a 
magnetron characterized by a narrow mode width and critical 
dependence on the voltage rate of rise of the modulating pulse. 

A hard tube modulator with rapid ri '■ capability and stable 
output was therefore selectee* magnetron. The 

incut for this modulator stage is obtained from a series of 
avalanche transistors which are driven into their avalanche 
mode of operation by the trigger input from the PRF generator. 
The rate of rise of the pulse output of the avalanche circuitry 
is controllable to some degree by adjustment of the terminating 
impedance. The pulse width is controlled by the length of the 
delay line whose charge feeds the energy into the avalanche 
transistors. The avalanche output of roughly 150 volts exhibits 
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typical rise times of 2 nanoseconds. Modulator pulse widths* 
from 20 to 60 nanoseconds have been obtained. 


The hard tube modulator approach U3ing the avalanche 
transistor driver has resulted in stable and repeatable RF 
operating conditions, requiring little or no modulator voltage 
adjustments on a day-to-day basis. The additional volume and 
power required by the avalanche circuitry is relatively in¬ 
significant in comparison to the improvement in operating 
reliability. 

System amplitude stability was verified by boxcarring 
the video return from a solid, isolated target (tall smoke 
stack) located at 31 microseconds in range, see Figure 44. 

A low frequency analysis performed on this boxcarred signal 
output failed to show any low frequency scintillation content 
which could be attributed to instabilities in the radar. A 
similar low frequency spectrum analysis was carried out on the 
returns from bare trees, see Figure 39, at roughly 2 micro¬ 
seconds in range. It disclosed a low frequency scintillation 
content extending to only 25 Hz. This is considerably less 
than the 300 Hz extent of scintillation frequencies obtained 
in spectrum analyses performed against trees in full foliage 
at Ballistics Research Laboratories, Aberdeen, Maryland in 
the summer of 1965. (See ECOM-01253-2, Appendix C) 


1.4.2 Experimental Results 

H series of preliminary experiments have been conducted 
using the experimental radar described in section 1.4.1. The 
data was obtained at the nominal pulse width of 60 nanoseconds 
and at the narrower pulse width of 20 nanoseconds. 

Figure 37 shows the 60 nanosecond transmit*-' pulse at 
a power level of 900 watts (peak). The shape < f. pulse is 
quite square. The ringing is due primarily to ...» measurement 
circuitry. 

Figure 38 shows the return from a large tree with no 
foliage. The turn is quite strong and re^aui.vely steady. 

Figure jho :he power spectral density of thir return. Note 
also the steady peaked return of Figure 38 which subtends approxi¬ 
mately 200 ns., i.e., 100 feet, in range. This is much larger 
than the tree dimension. 

Figure 40 shows the discrete return from an aluminum 

cable. 


Figure 41 shows the return from a large 
has the distinguishing twin peaks as in Figure 


tree, 
38. 


and again. 
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Figure 37 


60 ns Pulse, Transmitted 
Power of 900 watts (Time 


at a Peak 
scale * 20 ns/cm) 



Figure 38. Return from a Tree with No Foliage 
(Time scale = 200 ns/cm) 
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Figure 39. 


Clutter Spectrum from Trees at 1000 Feet 
in Range 
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Figure 40. Return from an Aluminum Cable 
(Time scale = 200 ns/cm 



Figure 41. Return from a Tree with No Foliage 
(Time scale = 200 ns/cm) 

















Figure 42 shows the return from a wooden house, or 
more explicitly, from the metal gutters on the house. 

Figure 43 shows the return from a 30-foot telephone pole. 
The return is extremely steady. A similar return has been 
observed from a brick chimmney at a range of 2.9 miles. The 
width of the return in each case is essentially the pulse width, 
60 ns. The power spectrum of the chimney is shown in Figure 44. 
The soectrum is almost white with very little low frequency 
scint; ' ’.at ion. 

'•'iuurc 45 shows the narrow pulse shape, 20 nanoseconds. 

Figure 46 shows the return from a house using the 
narrow pulse. Note the marked decrease in clutter level com¬ 
pared to Figure 42. 

A series of experiments, to obtain data from a variety 
of cables and towers, is planned for the next quarter. 


1.4.3 Probability Densities of Experimental Data 
1.4.3.1 Introduction 


A program of millimeter-wave radar experiments, initiating 
our investigation ot the radar returns from clutter targets, was 
carried out at the Ballistics Research Laboratories, and is fully 
described in Appendix C of Technical Report ECOM-01253-2, March, 
1966. 


We are presently concerned with obtaining a description 
of such video returns in terms of their power spectral density 
and probability density function. The boxcarred, radar-video 
signal under investigation was either recorded with a Precision 
Instruments tape recorder for later analysis with respect to 
the signal's probability density function or analyzed immediately 
with respect to its power spectral density by means of a Panoramic 
Subsonic Analyzer. The results of the power spectral analysis 
were renorted in the previously referenced technical report. 

Our purpose here, is to study the probability density function 
of the magnitude of the clutter and composite radar returns. 
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Figure 42. Return from the Metal Gutters 

on a Wooden Hou^e (Time scale = 200 ns/cm) 



Figure 43. Return from a Telephone Pole 
(Time scale = 200 ns/cm) 














Figure 44. Radar Spectrum of Return from Chimney at 
2.9 Miles in Range 
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Figure 45. 20 ns Pulse (Time scale - 20 ns/cm) 



Figure 4€. Return from the Metal Gutters on a 

Wooden House (Time scale = 100 ns/cm) 


















1.4.3.2 Data Processing Techniques 


The data processing work required to characterize the 
video signal in terms of the probability density function of 
its magnitude was performed at the United Aircraft Corporation 
computer facility in East Hartford, Connecticut. A block 
diagram of the system used to obtain an estimate of the probabil¬ 
ity density function is shown in Figure 47. The previously 
recorded, boxcarred video signal was sampled by the computer at 
a 1 kHz rate and sorted according to magnitude, into 200 bins, 
each 20 millivolts in width. The number of samples in each bin 
was then plotted out by the X-Y recorder. A typewriter printout 
of tne data was also obtained. The entire operation proceeded 
under the supervision of the control unit. 

A test run was performed in which more than 10** samples 
were taken of the signal 

y = f (= sin (6) (27) 

with 6 = t and u = 1, at a 1 kHz rate. The experimentally 
determined density function, p(y), is compared, in Figure 48, 
with the theoretically predictable density function 

p 2 (y) = 1 — [13] (28) 

►1-y 2 

The area under the experimentally determined density function 
has been normalized to unity. The two densities are seen to 
be in good agreement. 


1.4.3.3 Density Measurements of Radar Video 

Figures 49 through 54 are reproductions of the outputs 
of the X-Y recorder for six different radar signals being 
studied. The numbers in parenthesis refer to figure numbers 
from Appendix C of Technical Report ECOM 0-1253-2 in which the 
corresponding power spectral density is displayed. Figures 49 
and 50 are the probability density functions of clutter only, 
at two different wind velocities. Figures 51 through 53 are 
the probability density functions of isolated, stationary targets 
and Figure 54 is the probability density function of a composite 
target. 


The area under each of the experimental curves in Figures 
49 through 54 is the same, since the total number of samples 
taken in each case was 2 x 10 1 *. In order for the figures to 
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Figure 47. Block Diagram of the System Used to 
Obtain Estimates of the Probability 
Density Function of the Magnitude of 
a Radar Echo 
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Figure 49. Probability Density Function of Radar Echo 
from a Line of Trees Beginning at 560 Yards 
Range with 0-3 mph Wind Velocity 
(Ref. ECOM - 01253-2, figure C-4) 
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Figure 50. Probability Density Function of Radar Echo 
from a Line of Trees Beginning at 560 Yards 
Range with 10-14 mph Wind Velocity 
(Ref. ECOM - 01253-2, figure C-5) 
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Figure 51. Probability Density Function of Radar Echo 
from Reference Sphere 
(Ref. ECOM - 01253-2, figure C-6) 
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Figure 52. Probability Density Function of Radar Echo 
from Isolated Telephone Pole 
(Ref. ECOM - 01253-2, figure C-8) 
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Figure 53. Probability Density Function of Radar Echo 
from Vertically Suspended Cable 
(Ref. ECOM - 01253-2, figure C-7) 
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properly represent probability density functions, this area has 
been normalized to unity. In each case, a Gaussian probability 
density function 

g 2 (y) = ——r - exp [- | (-^-) 2 ] (29) 

0 /2 IT 


has been drawn on the same set of axes as the experimentally 
determined density function, gj (y). The Gaussian density is 
completely specified by specifying its mean, m, and standard 
deviation, o. The mean of the Gaussian process was set equal 
to the computer calculated value of the sample mean of ^,ne 
video signal. In presenting the data the sample mean whs set 
equal to zero, in each case. The standard deviation was chosen 
in such a way that the area under the experimental curve from 
(m - o) to (m + c) is equal to the area under g 2 (y)/over the 
same range. ■ 


1.4.3.4 Conclusion 


Based on Figures 49 through 54 several tentative con¬ 
clusions can be made concerning the relationships between the 
nature of the radar target and the probability density function 
of its boxcarred, radar-video return. These are presented as 
follows: 

a. Each of the experimentally determined probability 
density functions can be well approximated by a 
normal distribution with a standard deviation 
selected as indicated in section 1.4.3.3. This 
result provides additional justification for our 
use of a normal approximation in the ground clut¬ 
ter simulator discussed in Appendix G of Technical 
Report ECOM-01253-2. 

b. A certain asymmetry about the sample mean can be 
observed to exist for each of the experimentally 
derived density functions. The asymmetry is more 
pronounced for clutter targets than for composite 
targets and for composite targets more than it is 
for isolated, stationary targets, with the exception 
of the reference sphere. Theoretical calculations 
based on a random scatter model 11] indicate the 
existence of such a trend. For low fixed-to-fluc- 
tuating target ratios, as in the case of clutter 
targets, the predicted distribution is an asymmetric 
Rayleigh distribute n. Whereas, for high fixed-to- 
fluctuating target ratios, as in the case of iso¬ 
lated stationary targets, the predicted distribution 
is a symmetric Gaussian distribution. Hence, the 

.ncrease in the symmetry of the density function 
as the fixed-to-fluctuating target ration increases. 
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c. Figure 49 is the probability density function of 
the return from a patch of trees moving in a 0-3 mph 
wind. Figure 50 is the probability density function 
of the same patch of trees moving in a 10-14 mph 
wind. The standard deviation of the Gaussian den¬ 
sity of Figure 49 was found to be equal to 0.58. 
Whereas, the standard deviation of the Gaussian 
density of Figure 50 was found to be equal to 0.40. 
Based on these two densities, it would appear that 
the higher the wind velocity, the narrower the 
probability density function. 

This result is difficult to explain on the basis of 
a random scatterer model. In such a model, it is 
assumed that the radar echo is the vector sum of 
the radar returns from each individual scatterer. 

If it is assumed that the magnitude of the radar 
return from each individual scatterer is an inde¬ 
pendent random variable, obeying a particular 
density law, then the probability density function 
of the overall return is a complex function of the 
magnitudes of the individual returns and the density 
functions which they obey. 

d. Figure 54 is the probability density function of 
a composite target consisting of the patch of 
trees of Figure 49 and the telephone pole of Fig¬ 
ure 52. The inclusion of a stationary target in 
the illuminated patch of clutter is seen to in¬ 
crease the probability of obtaining a sample of 
the video return whose magnitude lies in the 
neighborhood of the mean. 

It is important to note that the statements enumerated 
above are made on the basis of the limited number of pieces of 
experimental data presented here. The tentative conclusions 
presented can, however, serve to indicate areas of investigation 
for a series of experiments presently being conducted by means 
of the radar system described in section 1.4.1 of this report. 
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